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The Yang-Mills Schrodinger equation is variationally solved in Coulomb gauge for the vacuum 
sector using a trial wave functional, which is strongly peaked at the Gribov horizon. We find the 
absence of gluons in the infrared and also a confining quark potential. 
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1. Introduction. One of the most challenging prob- 
lem in particle physics is to explain the confinement of 
quarks and gluons in QCD. Several confinement mecha- 
nisms have been proposed in the past, the most promi- 
nent of which are perhaps the condensation of magnetic 
monopoles (dual Meissner effect) and center vortex con- 
densation (for a recent review see ref. 1]). Evidence for 
both mechanisms has been found in lattice calculations. 
However, the center vortex picture is perhaps more ad- 
vantageous in the sense, that these vortices can, in prin- 
ciple, be defined in a gauge invariant way [2| , while mag- 
netic monopoles arise as gauge artifacts after Abelian 
projection and are merely manisfestations of topological 
defects of the underlying gauge fields Q. Yet another 
confinement mechanism was proposed by Gribov [J] and 
further elaborated in ref. 5j. This mechanism, which is 
based on the infrared dominance of the field configura- 
tions near the Gribov horizon in Coulomb gauge, is com- 
patible with the magnetic monopole and center vortex 
pictures of confinement, given the fact, that these field 
configurations lie on the Gribov horizon [6(. 
In this letter we will explore Gribov's confinement mecha- 
nism by studying the vacuum sector of SU(2) Yang-Mills 
theory in Coulomb gauge in the Hamilton approach us- 
ing the variational principle. With an appropriate, phys- 
ically motivated ansatz for the Yang-Mills wave func- 
tional, which accounts for the dominance of the field con- 



figurations on the Gribov horizon, we will be able to de- 
scribe simultanously quark and gluon confinement. 
In previous studies of the Yang-Mills Schrodinger equa- 
tion in Coulomb gauge 0, M, , a different trial wave func- 
tion was used and also the non-trivial metric of the orbit 
space, induced by the Faddeev-Popov determinant was 
not fully included. In ref. 7] the curvature of orbit space 
was completely neglected, while in ref.0 its contribution 
to the gluon self-energy was partially omitted in the gap 
equation. We will find, however, that the proper inclu- 
sion of the Faddeev-Popov determinant is absolutely nec- 
essary to produce simultaneously quark and gluon con- 
finement. 



2. Yang-Mills theory in Coulomb gauge. The Hamilton 
approach to gauge theory is based on the Weyl gauge 
Aq = 0, in which the dynamical degrees of freedom are 
the spatial components of the gauge field, A(x). This 
gauge leaves still invariance under spatial gauge trans- 
formations. The latter can be fixed by implementing the 
Coulomb gauge dA = 0. In the Coulomb gauge the un- 
physical longitudinal degrees of freedom of A(x) can be 
completely eliminated by explicitly resolving Gaufi' law 
Z)n|\I r ) = Pext\^)i resulting in the following Hamilto- 
nian of the physical transversal degrees of freedom A^~ 
(diAt = 0) 
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Here g is the Yang-Mills coupling constant, II" (x) — 
i5/SAj-(x) is the canonical momentum operator conju- 
gate to the transversal gauge field A 1 - (x) , D = d + A 1 - 
is the covariant derivative in the adjoint representa- 
tion (A ab = f°-cb A c jabc being the structure constant) 



and J = Det(-dD) is the Faddeev-Popov determi- 
nant. Given the fact, that the Faddeev-Popov kernel 

(—dD) represents the metric tensor in the color space 
of transversal gauge connections Af- the first term in 
the Hamiltonian is the corresponding Laplace-Beltrami 



operator and gives the electric part of the Hamilto- 
nian. The second term gives the magnetic energy with 
Bj^l" 1 "] = ^€ijk[Dj, Dk] being the color magnetic field. 
This term represents a potential for the transversal gauge 
field. Finally the last term is the so-called Coulomb term, 
where p a = Af- T[\ +p% xt is the non-Abelian color charge 
density, which contains besides the gluonic part also a 
contribution from external quarks p ex t- 
In this letter we solve the Yang-Mills Schrodinger equa- 
tion Hty = E^S for the vacuum sector by the variational 
principle using the following ansatz for the Yang-Mills 
vacuum wave functional 
TV 



*[A ] = 



• J d 3 xd 3 x' A ± 1(x)uj(x,x')A ±a i (x') 
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where TV is a normalization constant and the kernel 
u(x,x') is determined by minimizing the energy. The 
wave functional is strongly peaked at the Gribov hori- 
zon, where the Faddeev-Popov determinant vanishes and 
thus reflects the fact, that the dominant infrared configu- 
rations, like center vortices or magnetic monopoles, lie on 
the Gribov horizon |6(. Furthermore, the wave functional 
being divergent on the Gribov horizon, identifies all 
configurations on the Gribov horizon, in particular those 
which are gauge copies of the same orbit. This identi- 
fication is absolutely necessary to preserve gauge invari- 
ance. In addition it topologically compactifies the (first) 
Gribov region. Thus the pre-exponential factor [A] 
drastically changes the properties of the vacuum state 
compared to those of a pure Gaussian, which was used 
in refs.Q, Q- In the case of QED, where the Faddeev- 
Popov determinant becomes a constant eq. (0 represents 

the exact vacuum wave functional with w(fc) = \/k 2 be- 
ing the energy of a free photon with 3-momentum k. 

3. S chiving er- Dyson equations. In the eval- 
uation of the vacuum energy (^\H\^) = 
J DA ± J{A ± }^*lA ± ]H-$i[A ± } the following ingredi- 
ents are required: 

1. The ghost propagator in the vacuum, which is defined 
by 

G=(n-d£>)- 1 \y) = (-3 2 )- 1 - , (3) 

9 

where d denotes the ghost form factor, which measures 
the deviations of the ghost from a free massless field. 
Evaluating this expectation value in the so-called rain- 
bow ladder approximation, where the self-energy of the 
ghost is given by the diagram shown in fig^, one obtains 
the following integral equation for the ghost form factor 
d in momentum space 



1 



d(k) 

where (Nc number of colors) 
d 3 q 



--ld(k) 
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FIG. 1: (a) The ghost self-energy in the rainbow ladder ap- 
proximation. The full (wavy) line represents the full ghost 
(gluon) propagator, (b) Ghost loop contribution to the gluon 
self-energy. 
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FIG. 2: Diagrammatic illustration of the integral equation 
for the Coulomb form factor / represented by the fat dot with 
two entries for external legs. 



2. The Coulomb form factor / defined by 

(y\(-d5)- 1 (-A)(-d5)- 1 \v) := g(-a)/g . (6) 

The calculation of this form factor consistent with the 
evaluation of the ghost form factor yields the following 
integral equation 



f(k) = l + I f (k) 



Nr 



= / d*q[l-(kq) 



d 2 (k-q)f(k-q) 
(k - q) 2 uj(q) 



(7) 



(8) 



which is diagrammatically illustrated in fig. 2. In this 
equation we will replace the full ghost form factor by its 
bare value d(k) = 1, in order to carry out the calculations 
consistently to 1-loop order. 

3. The scalar "curvature" of color orbit space 



1 1 
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m , (9) 



which is the ghost loop part of the gluon self-energy and 
represents the ghost part of the color dielectric suscep- 
tibility of the Yang-Mills vacuum. To 1-loop order this 
curvature is given by the diagram shown in fig. lb and 
equals 



X(k) = I x (k) , 
Nc 
4 ./ (2tt)3 



d(k - q)d(q) 
(k - q) 2 



(10) 



Evaluating the expectation value of the Hamiltonian JTJ 
with the trial wave functional to 1-loop order varia- 
tion of the energy leads to the following gap equation for 
the kernel u)(x, x') in the wave function J5J 



u 2 {k) 



X 2 (k) + lL 2 Hk) + 2x(k)li 1 Hk)+C 



(11) 
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where I® is an irrelevant constant, which drops out after 
renormalization, and 



d{k-qff{k-q) 



Nc 

4 J (2tt)3 r 1 v "» j (k - q y 

(u>(q)- X (q)) n -(u>(k)-x(k)y 



uj(q) 



(12) 



Equations Q), JJJ), l(10|l and 1)11(1 represent four coupled 
Schwinger-Dyson type of equations for the ghost form 
factor d(k), the Coulomb form factor f(k), the curva- 
ture x(^) an d the gluon energy uj{k). These equations 
contain divergent integrals and require thus regulariza- 
tion and renormalization. Fortunately, the asymptotic 
infrared and ultraviolet behaviour of the solutions does 
not depend on the details of the renormalization proce- 
dure used. 

4- Asymptotic behaviour. One can solve the coupled 
Schwinger-Dyson equations analytically in the ultraviolet 
k — > oo in the so-called angular approximation 10] . One 
finds then the following asymptotic ultraviolet behaviour 



d(k) 



X(k) 

Lj(k) 

f(k)- 



1 



i 



(13) 
(14) 



where \i is an arbitrary parameter of dimension mass. 
The first equation means, that the gluons behave asymp- 
totically like free particles, while the second one implies 
that the space of gauge connections becomes asymptoti- 
cally flat. The ghost and Coulomb form factors, d(k) and 
f(k), deviate asymptotically from that of a free massless 
field by the anomalous dimension factor 1 / y/lnk/fx. 
In the infrared one can rigorously show, that x(fc — » 0) = 
uj{k — > 0) [ill. In addition, adopting the angular approx- 
imation [l(j| one finds the asymptotic solution for k — > 



u(k) = x(k) 



1 

k 



d{k) 



, /(fc-0) = l. (15) 



The first relation implies, that in the infrared the gluon 
energy diverges and equals its self-energy part generated 

by the ghost loop, i.e. its free part "vk~^ has dropped out. 
The infrared diverging gluon energy is a manifestation of 
gluon confinement. It implies an infrared vanishing gluon 
propagator, which violates positivity, and accordingly the 
gluons do not occur as asymptotic particle states in S- 
matrix |12| .|14|. Furthermore, the infrared diverging 
ghost form factor and the infrared finite Coulomb form 
factor f(k — ► 0) = 1 is precisely the infrared behaviour 
needed to produce a linearly rising confining potential. 

5. Renormalization. To regularize and renormalizc 
the divergent Schwinger-Dyson equations we use a 3- 
momentum cut-off and a momentum substraction scheme 
similar to the one used in refs.0, @- However, in the 



present case new features arises in the renormalization 
of the gap equation ((TTl) due to the full inclusion of the 
curvature ©■ The details will be given elsewhere [Tlj . 
After renormalization one obtains the following set of 
Schwinger-Dyson equations 



m = W) AId{k) 

X (k) = *(/*) + AI x (fc) 
f(k) = f(ji) + M f (k) 



(16) 

(17) 
(18) 



co 2 (k) = k 2 -^ + (AI x (k)) 2 +CAI x (k) + Ali 2 \k) 
+2[ x ^) + A/ x (fc)]A/W(fc)+ w 2 ( M ), (19) 

where \i is the renormalization scale and we have intro- 
duced the abbreviations 



AJ d (fc)=]im (I d (k,A)- J d (/x,A)) etc 



(20) 



Furthermore u(fi), d(fi), x(m)> /(/«) an d £ = 2[x(/i) + 
Iu (a*)] are renormalization constants, which were de- 
termined as follows: u>(n) is used to fix the energy scale 
and drops out from the Schwinger-Dyson equations by 
rewriting the latter in terms of dimensionless quantities. 
d(n) enters only the Schwinger-Dyson equation (|16fl for 
the ghost form factor, which does not contain the re- 
maining renormalization constants. The ultraviolet be- 
haviour of d(k) found above in ea.l|14|) is independent 
of d(fi) but the infrared behaviour of d(k) depends cru- 
cially on d(ji). As long as d(ji) is smaller than some 
critical value d cr , d{k) approaches a (finite) constant for 
k — > 0. At the critical value d(fi) = d cr the ghost form 
factor d(k) diverges for k — > and above the critical value 
d(p,) > d cr no solution for d(k) exists. This critical value 
is the only value, which produces the infrared diverging 
ghost form factor found above analytically (|15l) . Fur- 
thermore in D = 3 (which will be considered elsewhere) 
a self-consistent solution to the coupled Schwinger-Dyson 
equations exists only for this critical value. Therefore we 
choose d(fi) = d cr . Fortunately, the self-consistent solu- 
tion is quite insensitive to the remaining renormalization 
constants x(aO an d £, which we have chosen for the defi- 
niteness as x(a*) = 0: fit 1 ) = 1 an d ^ = 0. 

6. Numerical Results. The renormalizcd coupled 
Schwinger-Dyson equations eg . ifTBT) . l(T7jl . (|T%|) and (fT^f) are 

solved by iteration (without resorting to the angular ap- 
proximation) and the results are shown in figs. 3 and 4. 

All these numerical results are in full accord with the 
ultraviolet and infrared behaviour extracted above ana- 
lytically within the angular approximation. Finally, fig. 5 
shows the static quark potential which is obtained in the 
present approach as the expectation value of the last 
term in the Yang-Mills Hamiltionian, cq.Q, when the 
color density p a is identified with that of a static quark- 
antiquark pair. The obtained potential interpolates be- 
tween a Coulomb potential at small distances and an (al- 
most) linearly rising confinement potential at large dis- 
tances. Numerically we find, that its Fourier transform 
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FIG. 3: The ghost form factor d(k) (full line) and the 
Coulomb form factor f(k) (dashed line). 
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FIG. 4: The gluon energy cj(fc) (full line) and the curvature 
X(k) (dashed line). 

diverges for k — > as 1/k 3 instead of 1/fc 4 , which is 
required for a strictly linear rising potential |l5j . Ap- 
proximating the confining potential by a linear one we 
can fix the energy scale by the string tension. In ref.0 it 
was found that the string tension in the Coulomb poten- 



tial, cr cou i, is by about a factor of three larger than the 
asymptotic string tension a extracted from large Wil- 
son loops. Using a cou i = 3cr and the canonical value 
a = (AAQMeV) 2 one finds that the minimum in the gluon 
energy (see fig. 4) occurs at k m i n « lAGeV, which is of 
the order of the glue ball mass. This corresponds to a 
minimal single gluon energy u>(k m i n ) w 3GeV. 
To summarize, by approximately solving the Yang-Mills 
Schrodinger equation in Coulomb gauge by means of the 
variational principle using the trial wave functional J2J, 
which embodies the dominance of the field configurations 
on the Gribov horizon, we have been able to describe 
simultaneously quark and gluon confinement. A more 
detailed presentation will be given elsewhere. 
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FIG. 5: The static quark potential. 
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